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A SINGLE-COLLISION MODEL FOR ELECTRON-BEAM 


CURRENTS BETWEEN PLANE ELECTRODES* 
by Charles M. Goldstein and Arthur W. Goldstein 
Lewis Research Center 

SUMMARY 

Space-charge and potential distributions are calculated for the region between two 
plane electrodes. The effect of elastic, hard- sphere collisions between the electrons of 
the injected beam and the neutral gas molecules is incorporated into the analysis with the 
assumption that the electrons will suffer at most one scattering collision. The injected 
beam, itself, is assumed monoenergetic. Calculations of current as a function of voltage 
are compared with results from a multiple- collision model; very good agreement is found 
for an electrode- spacing to mean-free-path ratio of 0. 1, and a maximum deviation of only 
1. 5 percent was observed for a ratio of 0. 5. Current-voltage characteristics are also 
calculated for various values of emission current and electrode spacing. Collisions are 
found to have a pronounced influence on potential and density distributions, current, and 
region of time- independent operation. 


INTRODUCTION 

An important problem in the study of gaseous electronics and ionic phenomena is to 
find a means of analyzing the effect of collisions for those circumstances wherein they 
cannot be neglected but where a diffusion approximation is not yet applicable. This inter- 
mediate regime is of major importance in the low density plasma diode, the plasma 
sheath, the ion rockets, and the cross-sectional measurements. Presented herein is a 
simple, theoretical first approximation to the effect of collisions in this regime. 

The particular problem treated is that of electron flow between plane parallel elec- 
trodes (emitter and collector) in the presence of an external electric field (see fig. 1) 
and collisions with a neutral target gas. The emitted electron beam is assumed mono- 
energetic and directed normal to the electrode surfaces. In addition, the following sim- 
plifying assumptions are made: the neutral gas particles are of infinite mass (perfect 

The material contained herein was presented at the Seventh International Conference 
on Phenomena in Ionized Gases, Belgrade, Yugoslavia, August 22-27, 1965. 
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Lorentzian gas); electron- neutral collisions 
are elastic, hard- sphere collisions; and, the 
electrons suffer at most one scattering colli- 
sion. The study is further restricted to the 
regime where a time- independent solution can 
be considered valid. 

This single- collision model should provide 
good qualitative characteristics for a dilute 
gas; it will also provide a much-needed check- 
point for all future exact or approximate solu- 
tions of the more accurate physical models. 
The pursuance of this problem was in no small 
way influenced by the need for just such a 
checkpoint for a Monte Carlo problem (ref. 1) . 
Not only did this model provide the needed checks for the Monte Carlo solution, but these 
latter solutions demonstrate, in turn, the range of validity of the one- collision approxi- 
mation. A comparison of the results obtained from the two models will be discussed. 

A review of the electron vacuum diode phenomena, which are of interest in the pres- 
ent study, is given in the following section. A more extensive review of research on 
monoenergetic electron flow in a vacuum diode is to be found in reference 2. 


Figure 1. - Configuration of one-dimensional field, flux, and 
electrodes. 


REVIEW OF COLLISIONLESS RESULTS 

Figure 2 is a typical current-voltage characteristic obtained from a time- independent 
analysis (ref. 3) for the case of no collisions; figure 3 shows the corresponding potential 
distributions. As the collector voltage is decreased from A in figure 2, the full emis- 
sion current reaches the collector until point B is reached. The curve labeled B in 
figure 3 represents the potential distribution corresponding to point B in figure 2. If 
the collector potential is now decreased infinitesimally, the potential minimum drops 
discontinuously to zero as indicated by curve B’ in figure 3. This situation corresponds 
to point B' on the current- voltage characteristic (fig. 2). The point where the minimum 
potential and its slope are both zero is called a "virtual cathode. " At this point, the 
ideal monoenergetic beam has zero velocity (infinite charge density). Other values of 
collector potential are determined by the amounts of the incident electron current that 
are reflected or transmitted at the virtual cathode. If the collector potential is now in- 
creased from, say, point C (fig. 2), a virtual cathode will continue to exist until point D 
is reached where the potential undergoes a discontinuous change and a normal potential 
minimum is once again formed (curve D', fig. 2). 
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Figure 2. - Current-voltage characteristic; no collisions. 



Figure 3. - Potential distributions; no collisions. 


The hysteresis predicted by the collisionless current-voltage characteristics were 
not observed experimentally (ref. 4). More recently, the problem was simulated on a 
high-speed digital computor (ref. 5). The results indicate that a virtual cathode is never 
formed. Instead, for anode potentials below the critical potential (point B, fig. 2 and 
curve B, fig. 3), the potential minimum appears to oscillate nonuniformly with a period 
of roughly the order of the electron plasma frequency corresponding to the electron den- 
sity at the potential minimum. Thus, in effect, only time-dependent solutions appear to 
exist below the critical potential. 


ANALYSIS 

Solutions to the single- collision problem are obtained by iteration in the following 
manner: an approximate potential distribution is first assumed for a given collector po- 
tential; the contribution of the electrons in the primary beam to the space charge is then 
easily calculated; the contribution of the scattered electrons is obtained from a solution 
of the Boltzmann transport equation; the total space charge is then employed, with 
Poisson’s equation, to obtain a more accurate potential distribution. The entire process 
is iterated until sufficient accuracy is obtained. 

Poisson’s equation for the interelectrode region in dimensionless variables becomes 

<p”(y) ~ Cn(y) (1) 
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(All symbols are defined in the appendix. ) The parameter C is proportional to the 
strength of the space-charge effects. A low value of C indicates weak space-charge ef- 
fects and vice versa. 

The dimensionless charge density n(y) may be expressed as 


n(y) = n p (y) + n g (y) (3) 

where n p (y) and n g (y) are the contributions to the charge density at y from the primary 
and scattered electrons, respectively. 


Density of Electrons in Primary Beam 

The assumption of hard- sphere collisions and a constant mean-free-path X implies 
(see ref. 6) 

N p(x)v p (x) = ^ e' x / A (4) 

where N (x) is the density and v (x) is the velocity of electrons in the primary beam 

lr r 

(i. e. , those electrons that have not been scattered). In this equation 
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( 5 ) 


v p <x) = 


and, therefore, from equations (2), (4), and (5), 


n p (y)= e 


■KY 


Vi + <p(y) 


( 6 ) 


where 


X 


( 7 ) 


Equation (6) represents the density distribution of electrons in the primary beam. 


Density of Scattered Electrons 

Source intensity . - To solve Boltzmann’s equation for the scattered electrons, the 
source intensity of scattered electrons is equated to the rate of loss of electrons from 
the primary beam. Equation (4) represents the flux of the primary beam electrons as a 
function of distance between the electrodes. The rate at which electrons are lost from 
the primary beam in distance dx is obtained by differentiating equation (4) with respect 
to x in the following way: 


t N p (x> V x) ' - 


N p( x )v p (x) 

X 


Hence the rate of loss of electrons from the primary beam per unit volume in phase 
space (dr dv) is 


N (x)v (x) 

P ~ P — 8 (v x - Vp)5 (v y ) 6 (v z ) (8) 

Since hard- sphere scattering is energy conserving and isotropic, the rate of gain of 
scattered electrons per unit volume of phase space can be written in the form 
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(9) 


R(x,v)s(v 2 - v 2 ) 


When the gains and losses are equated (by integrating the integrals of eqs. (8) and (9) over 
the entire velocity space), there results 


N„(x) 

R(x, v ) = R(x) =-^ — 

2ttX 


( 10 ) 


Expressions (9) and (10) define the source intensity of the scattered electrons in phase 
space. 

Equations for scat tered electrons. - The Boltzmann transport equation for the scat- 
tered electrons with the derived collision term is 

V I » “V (11) 

x 3x m dx dv x 2ir\ ' 

where f(x, v ) is the velocity distribution function of the scattered electrons. Integrating 
equation (11) with respect to v y and v z and nondimensionalizing result in 

u ag < y > u) + A <p'( y) ASM = s(y) (12) 

3y 2 3u 


where 
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(13) 
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The dimensionless density of scattered electrons n (y) is given by 

f n max 

n g (y) = J g(Y,u)du (14) 

u min 

where g(y, u) is a solution of equation (12) in the range bounded by u m ^ n and u max 
(which are themselves functions of y) . 

Solution of Boltzmann equation . - Equation (12) will be solved by the method of char- 
acteristics (ref. 7) . After a collision has occurred, the normal velocity components Vy 
and v_ of each scattered electron are constant (no further collisions?); hence, the quan- 
tity 

E = u 2 -</>(y) (15) 

becomes a constant of motion of the scattered electrons. Transforming from the (y, u)- 
space to the (y, E)- space where 


u = ± Ve + cpi y) 


(16) 


results in equation (12) becoming 


u(y, E) iS (y, E) = S(y) (17) 

3y 

Integration of equation (17) along constant E from a lower limit ^ to y gives 


g(y, E) = g(£ , E) 



dt 


(18) 


The bookkeeping becomes easier if the intermediate distribution functions gj(y, E) and 
g 2 (y, E) are employed 


gj (y, E) = g(y < y m , E) = gj(y, E) + gj(y, E) 
g 2 (y, E) = g(y > y m , E) = g£(y, E) + g 2 (y, E) 


(19) 


where y m is the location of the potential minimum, the g£(y, E) and g^(y, E) are the 
distribution functions for electrons with velocity u > 0 and u < 0, respectively. From 
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equations (18) and (19) 


gj(y, E) = g*£ , E) ± fj K(t, E)dt 


( 20 ) 


where 


K(t,E) = 


S(t) 


Ve + <p(t) 


The boundary conditions at the electrodes, if perfect absorption is assumed, are 


( 21 ) 


gj(0, E) = 0 (no scattering at emitter) j 
gg(l, E) = 0 (no scattering at collector)! 


( 22 ) 


Continuity of the distribution functions at the potential minimum is also demanded as 
shown in the following equation: 


Sluin’ = E ) 


(23) 


u 



(a) y < y m ; region 1. 



(b) y>y m ; region 2. 

Figure 4. - Phase space of scattered electrons. 


The integrals of equation (20) are to be evalu- 
ated along lines of constant E. 

The limits of integration may be readily ob- 
tained with the aid of figure 4, which represents 
the phase space of scattered particles with 
E, u-axes. Contours of constant position (y = 
constant) are parabolas given by (cf. eq. (15)) 

2 

u = E + cp( y) 

There are two regions in phase space; the first 
(fig. 4(a)) describes the region in physical space 
0 < y < y m where all the electrons are deceler- 
ated (in the vector sense), and the second (fig. 
4(b)) describes the region in physical space 
y m < y < 1 where all electrons are acceler- 
ated. In these diagrams, the paths of individual 
electrons are vertical lines (constant E) going 
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downward in region 1 (fig. 4(a)) and upward in region 2 (fig. 4(b)). The integrals of equa- 
tion (20) are taken along the same paths. 

The point of origin of a particle is determined by the location y of scattering and the 
axial velocity u after scattering. These determine the constant of motion E, which is 
always less than 1.0. The boundary conditions (eqs. (22)) imply g = 0 on y = 0, u > 0 
(fig. 4(a)) and on y = 1, u < 0 (fig. 4(b)). 

As an example, the trajectories of an electron starting near the emitter (represented 
by arrow 1 in fig. 4) shall be followed. This electron has a positive u-component of ve- 
locity as shown. Following a path of constant E toward the potential minimum, u de- 
creases. After passing the potential minimum, the electron is accelerated toward the 
collector (arrow 2 in fig. 4(b)) . Arrow 3 (fig. 4(b)) represents an electron with the same 
value of E but with a negative value of u resulting from a backscattering collision. 

The magnitude of u decreases until this electron traverses the potential minimum and is 
then accelerated into the emitter (arrow 4, fig. 4(a)). 

Figure 4(a) also indicates that all electrons having positive u and E < ~<P m in the 
region y < y m will be reflected by the potential field and returned to the emitter. Simi- 
larly, electrons with negative u and E < -<p m in the region y > y m (fig. 4(b)) will be 
reflected back to the collector. The trajectory of such an electron is represented by ar- 
row 5 in figure 4(b) . The value of y at which u = 0 is also indicated as y*. Physical- 
ly, this represents the turning point in the trajectory for a constant value of E due to 
the potential field. 

With equations (20), (22), and (23), the following may be written for the indicated re- 
gions of figure 4: 



s la< y ’ E > ■ 

sj(y, e > -<p m ) - ^ 

r y 

/ K(t, E)dt 
0 

(24) 

g lb (y, e) 

s g^(y,E > 

-?m> ' 82 < y m> E > 


+ f Ym K(t, E)dt 
y 

(25) 

gi c (y,E) 

= gj(y,E < 

"0 = h(y*,v< 


ry* 

+ / K(t, E)dt 

y 

(26) 


g ld (y ’ E) = 

gj(y, E < - Vm ) = ^ 

ry 

^ K(t, E)dt 

(27) 

g 2a (y> E ) 

s g 2 (y> E > 

-<?rr) = S^y m . E > 


r y 

+ / K(t, E)dt 

•'v 

(28) 


m 
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( 29 ) 


/* 1 

g 2b (y, E) = g 2 (y, E > -<p m ) = ^ K(t, E)dt 
g 2c (y, E ) s §2 (y ’ E < -O = j^ 1 K ^’ E)dt (so) 


g 2d (y> E) S g 2< y ’ E < “O = g 2< y *’ E < “O + 7" Y K(t > E)dt (31) 

y* 


Since the distribution functions must be continuous along constant E, 
be equated from equations (27) and (30), respectively: 

the following may 

gj(y* E < -<p m ) = gj(y*, E < -<p m ) = K(t, E)dt 

(32) 

g|( y *’ E < -<Pu) = §2 (y *’ E < -<Pn) = jf* K(tj E)dt 

(33) 

Hence, equations (26) and (31) become 


, v ry fj* 

gi c (y, E ) = J K(t, E)dt + 2 / K(t, E)dt 

1C * / 0 J y 

(34) 

g 2d (y, E) = 2 f y K(t, E)dt + f * K(t, E)dt 
*T* •'y 

(35) 


Substituting equation (33) in equation (25) and equation (32) in equation (28) results in 


g lb ( y , E) = / K(t, E)dt 


(36) 


g 2a (y, E) = fj K(t, E)dt 


(37) 


The limits on the integral in equation (14) can now be obtained in each region of the 
(y, E)- space by inspection of figure 4. For y < y , the respective contribution to the 
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densities becomes 


n la< y ) = 

/•Vl+<p(y) 

/ du y 

V<?(y)-<P m J 

f V 

f K [t, u 2 - (y)] dt 

0 

(38) 

n lb (y) = 

*>JUcp(y) 

/ dU J 

J yjcp(y)-cp m J 

/• 1 

f K [t, u 2 - <p(y)]dt 

y 

(39) 

fy/<p(y)-<p m *y 

n ic (y ) = / du / 

^0 •'0 

ry/<p(y)-<p m 

K[t,u 2 - <p(y)Jdt + 2 / du 

*'o ^ 

j-y* 

j K[t,u 2 -<p(y)]dt 

y 




(40) 

n i d (y) = J 

r ^/<p(y)-<p m 
i du i 

0 / 

r y 

1 K[t,u 2 - <p(y)]dt 
0 

(41) 

and for y > y m , 




n 2a (y) = 

. Vl-Kp(y) 

/ du 

\<p{y)-<p m J 

r y 

/ K [t, u 2 - <p(y)]dt 

0 

(42) 

n 2b< y > = 

/ .Vl+<^(y) 

/ du / 

J y/<p{y)-(p m A 

r l 

K [t, u 2 - <p(y)] dt 
T 

(43) 

” 2c (y) - 

f^<p(y)-<p m 

/ du 

J 0 

A 

/ K [t, u 2 - <p(y)]dt 

J y 

(44) 

/-V<p(y)-<p m 

”2d (y) = 2 / du / 

y 0 •'v 

.y 

K [t, u 2 - <p(y)] 

* 

/■V V (y)-?> m 

dt + / du / 

*'o A 

.1 

K[t,u 2 -<p(y)]dt 

r 


( 45 ) 
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A straightforward change of the limits of integration may be affected for equations 
(38), (39), (41), (42), (43), and (44), and for those integrals of equations (40) and (45), 
which do not contain y* as a limit. The integrals containing y* can be inverted, how- 
ever, by employing the following identity: 


f da f h(a,/3)d/3 = f d/3 f h(a,/3)do! 
*'0 * / 0 0 


(46) 


The second integral of equation (40) may be written in the form 


r (y) <Pm da /^hta,^ 
J 0 J 0 


(47) 


where 


yfa = u 

/3(t) a -<p(t) + <p( y) 
K[t, a - <p(y)] 


h (a,0) a 


Va /3'(t) 


j 


(48) 


(Note that it has been assumed that cp( t) is monotonic in the range of interest, hence 
/3*(t) and t may be expressed as functions of /3.) Application of the identity (eq. (46)) to 
equation (47) gives 


fVW-V m /■*&>-*’ m r y m r , 

/ d/3 / h(a, /3)da = 2 / dt / K[t, u - <p(y)J du (49) 

- n w /3 ** \7 \l rr, / tA _ OA 


><p(y)-<p 


V<p(y)- 


yj<p(y)-<p(t) 


The same procedure gives, for the first integral of equation (45), 

/•y r>/<p(y)-<p m 

2 / dt / K[t,u - <p(y)]dt (50) 

*y m yJ<p(y)-cp(t) 

After changing the order of integration in equations (38) to (45), the integration over 
u may be done in closed form. The indefinite integral is 
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(51) 



- <p(y)]du 



S(t)du 

- <p( y) + <p(t) 


= S(t)ln 


J 2 

u +fu - 


(p( y) + <p( t) 


Hence, integrating over u and collecting terms from equations (38) to (51) result in 


n s ( y ) 



dt • S(t) In 


V l + <?(y) + Vi + <?(t) 
V<p(y) - <? m W<?(t) - <? m 



dt • S(t) In 


VT + <?(y) +VT + <y(t) 
V^(y) - <? m ± V<p(t) - <p m 


(52) 


where the upper signs in the denominators pertain to y < y m and vice versa. 


Current to Collector 

To obtain the net collector current J, it is simpler to first obtain the net current to 
the emitter J^, 

— = f u-g 1 (0,E)du (53) 

J o J n<0 

or 

-p = y* <0 U [g lb (0, E) + g lc (0, E)] du (54) 

from equations (25) and (26). Substituting equations (34) and (36) into equation (54) and 
referring to figure 4 (p. 8) for the proper limits result in 
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Inverting the order of integration (using eq. (46) for first term in eq. (55)) and integrating 
with respect to u give 



-±<l+e-V 



(56) 


RESULTS 

The equations for the electron density distribution were solved simultaneously with 
Poisson's equation for the potential distribution and current. Solutions were obtained for 
several values of the space-charge parameter (C = 0. 1, 1, and 10) and the collision pa- 
rameter (L/A = 0, 0. 1, 0. 3, and 0. 5). For each pair of parameters, the collector poten- 
tial was varied from a large positive value down to the neighborhood of the critical poten- 
tial. 


Case 1 (C ■ 10) 

Figure 5 shows the overall current- voltage curves for high space- charge conditions 
and several values of L/A. For comparison, the results of the Monte Carlo calculation 
are shown, where multiple collisions are allowed but also limited to hard- sphere elastic 
collisions. The one-collision approximation yields collector currents that are accurate 
for L/A = 0. 1 and is accurate over most of the range for L/A = 0. 3, but it is somewhat 



Figure 5. - Current-voltage characteristics; C - 10. 
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Interelectrode potential minimum, <p m j n Interelectrode potential, <p(y) 



Interelectrode distance, y 


Figure 6. - Potential distributions near transition voltage; 
C = 10; L/X = 0. 5. 



^0 2.5 3.0 3.5 4.0 4.5 5.0 


Collector potential, ip{l) 

Figure 7. “ Critical potential and minimum potential with mean-free- 
path; C * 10. 


too high (by as much as 1. 5 percent) 
for L/X = 0. 5. These curves show 
that collisions cause a backscattering 
of electrons to the emitter even for 
large collector voltages (J/J Q remains 
less than unity) . 

Typical potential distributions 
(shown in fig. 6) are qualitatively the 
same as those that exist without colli- 
sions. The only apparent effect of col- 
lisions is to lower the potential mini- 
mum for a given collector potential. 

The effect of collisions on the 
critical collector potential and the cor- 
responding potential minimum is shown 
in figure 7. The critical potential is 
defined at the point where the slope of 
the curve of potential minimum plotted 
against collector potential becomes in- 
finite (ref. 4). It is observed that col- 
lisions reduce the range of stable oper- 
ation by increasing the critical value 
of the collector potential. In addition, 
the corresponding potential minimum 
decreases. 

The electron density distribution 
obtained by the one- collision method 
for C = 10 and L/X = 0. 5 is shown in 
figure 8. Curve C shows a distribu- 
tion corresponding to an electrical field 
of zero at the emitter (cp'( 0) = 0). In 
this case, the potential increases mono- 
tonically with y, the primary beam 
electrons accelerate, and their density 
decreases. Scattered electrons return- 
ing to the emitter increase the density 
there over the emitted value of 1.0. 
Density curve A corresponds to the 
lower potential distribution of figure 7; 
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the maximum density location occurs in the neighborhood of the minimum potential dis- 
tribution (viz. , the minimum u-component of electron velocity). The curve labeled B 
shows the density distribution (n^) of members of the primary beam for the same poten- 
tial curve (<p(l) = 2. 55). The density of scattered electrons (n g ) is seen to be in the 
range 0. 17 to 0. 29. 


Case 2 (C ■ 1) 

In this case the cur rent -voltage curves (fig. 9) are qualitatively similar to those for 
case 1. Here, however, the critical potential is much lower (<p( 1) « -0. 9). Note that as 
the collector potentials become very high, the saturation currents for this case approach 
the values for case 1 (fig. 5). As the collector potential increases, the effect of the 
space charge decreases; the only mechanism, then, that maintains the saturation current 
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Figure 10. - Potential distribution; C = 1; L/X = 0.5. 



Figure 11. - Density distribution; C s l; L/A = 0. 5. 

below the emission current (J/J Q = 1) at very high collector potentials is the scattering. 

The potential distributions for C = 1 and L/X = 0. 5 are shown in figure 10. As 
the collector potential cp( 1) is reduced in the sequence of curves D, C, B, and A, the po- 
tential minimum first approaches the collector, and then recedes, as indicated by the 
position of the circles in the figure. For this case the potential minimum never reaches 
the collector (i. e. , the slope of the potential never becomes zero at the collector). The 
corresponding density distributions are shown in figure 11, where it is observed that for 
this value of C the scattering begins to dominate the space charge as is evidenced by 
the lack of correlation in positions of potential minimum and density maximum. The 
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Collector potential, $o(l) 

Figure 12. -Current-voltage characteristics; C = 0. 1 

density of the primary beam is also 
shown for cp(l) = -0. 9 and for (p (1) = 0.7. 

Case 3 (C = 0. 1) 

Here, a qualitative change in the 
current- voltage curves (fig. 12) for col- 
lector potentials near cp( 1) = 0 and lower 
is observed. The collector potentials co- 
inciding with 9 ?’( 1 ) = 0 and <p‘(0) = 0 de- 
fine different regimes of operation as indi- 
cated by the corresponding potential dis- 
tributions, curves A and B, in figure 13. 
The region between these two curves is 
known as the space- charge region. As 
the collector potential is decreased from 
that of curve B to that of curve A, the po- 
tential minimum moves from the emitter 
to the collector. Below curve A the po- 
tential minimum is located at the collector, 
and above curve B the potential minimum 
is located at the emitter. The almost lin- 
ear potential distributions outside the 
space-charge region (curves C and D), 
and the small extent of the space- charge 



Interelectrode distance, y 

Figure 13. - Potential distribution; C = 0. 1; LA = 0. 5. 



0 .2 .4 .6 .8 1.0 

Interelectrode distance, y 

Figure 14. - Density distribution; 0 = 0. 1; LA = 0. 5. 
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region are indicative of weak space-charge effects (viz. , small values of C. ) 

The density distribution n(y) of electrons, corresponding to the potential distributions 
of figure 13, is shown in figure 14. Here the scattering and acceleration are easily shown 
to have opposite effects on electron density. The distribution A corresponds to a greater 
number of electrons being reflected back to the emitter than does distribution B. Hence, 
in general, the electron density of A is higher. Near the emitter, however, the return- 
ing electrons of curve A are accelerated to a higher mean velocity than those of curve B, 
due to the higher electric fields. The increased mean velocity more than compensates 
for the higher back current, so that a net decrease in density results. 


CONCLUSIONS 

The calculations show that adequately accurate prediction of the effect of collisions 
on diode current-voltage characteristics are obtainable from the single- collision model 
if the spacing L is sufficiently small compared with the mean free path X; errors of 
only 1. 5 percent in current were obtained with L/X = 0. 5. 

Backscattering of electrons causes a reduction in the current reaching the collector. 
This reduction becomes most pronounced at low collector potentials where the interaction 
between scattering and the space- charge effects are strongest. These effects are shown 
to be accentuated with increase in the magnitude of the space-charge parameter C. In 
addition, the scattering causes a slight reduction in the range of stable operation. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 14, 1965. 
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APPENDIX - SYMBOLS 


C space-charge parameter, 

eq. (2) 

E constant of motion defined in 

eq. (15) 

e electronic charge 

f(x, v ) velocity distribution function 

(v.d. f. ) of scattered elec- 
trons 

g(y, E) dimensionless v. d. f. for 

scattered electrons, eq. (13) 

h(a,0) variable defined in eq. (48) 

J net electron current to collector 

j net electron current to emitter 

ili 

J electron emission current 

o 

K(t, E) kernel defined in eq. (21) 

L interelectrode spacing 

m electron mass 

N(x) electron density 

n(y) dimensionless electron density 

defined in eq. (2) 

R(x, v) source kernel, eq. (10) 

S(y) dimensionless source intensity 

defined in eq. (12) 

u dimensionless x- component of 

velocity, eq. (13) 

V(x) potential with respect to the 

emitter 

v velocity, v x , v , v z 


coordinate normal to electrode 
surfaces 

dimensionless coordinate nor- 
mal to electrode surfaces 
defined in eq. (2) 

y* value of y where u = 0 in 

integration along constant E, 
shown in fig. 4 

a, /3 variable of integration, eq. (51) 

6(x) Dirac Delta- function 

k dimensionless mean-free-path 

defined in eq. (7) 

X mean-free-path (mfp) 

r volume 

cp(y) dimensionless potential defined 

in eq. (2) 

Subscripts: 

a, b, c, d regions of fig. 4 
m value at potential minimum 

max maximum 

min mimimum 

0 initial conditions (at emitter) 

p primary beam electrons 

s scattered electrons 

1 

2 y>y m 

Superscripts: 

+ flow toward collector 

flow toward emitter 


[Cgs units are used throughout. ] 
x 

y 
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